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ABSTRACT 


In weakly ferromagnetic materials, already small changes in the atomic configuration triggered by tem¬ 
perature or chemistry can alter the magnetic interactions responsible for the non-random atomic-spin 
orientation. Different magnetic states, in turn, can give rise to substantially different macroscopic prop¬ 
erties. A classical example is iron, which exhibits a great variety of properties as one gradually re¬ 
moves the magnetic long-range order by raising the temperature towards and beyond its Curie point 
of = 1043 K. Using first-principles theory, here we demonstrate that uniaxial tensile strain can also 
destabilize the magnetic order in iron and eventually lead to a ferromagnetic to paramagnetic transi¬ 
tion at temperatures far below t £. In consequence, the intrinsic strength of the ideal single-crystal 
body-centered cubic iron dramatically weakens above a critical temperature of ~ 500 K. The discovered 
strain-induced magneto-mechanical softening provides a plausible atomic-level mechanism behind the 
observed drop of the measured strength of Fe whiskers around 300-500 K. Alloying additions which 
have the capability to partially restore the magnetic order in the strained Fe lattice, push the critical 
temperature for the strength-softening scenario towards the magnetic transition temperature of the un¬ 
deformed lattice. This can result in a surprisingly large alloying-driven strengthening effect at high 
temperature as illustrated here in the case of Fe-Co alloy. 
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Iron is one of the most abundant elements in the Milky Way and constitutes the major component of 
earth’s core. 1 Processed iron is the main ingredient in steels and other ferrous materials, whose technol¬ 
ogy is a central pillar in today’s industrial world. Beyond the obvious practical interests, the physical 
properties of Fe have attracted great attention in many fields of sciences. Probably the most conspicuous 
phenomenon is the marvellous coupling between magnetism and mechanical properties. Although rec¬ 
ognized long time ago, 2 it is only recently that atomic-level tools can give insight into the mechanisms 
behind the magneto-mechanical effects in Fe and its alloys. Magnetic order in solid iron emerges with 
the formation of a crystalline lattice. The stable body-centered cubic (bee) phase (a-Fe) is characterized 
by long-range ferromagnetic (FM) order below T®. The face-centered cubic (fee) structure appears in 
the phase diagram at temperatures above 1189K in the paramagnetic (PM) state (y-Fc). At low tem¬ 
perature, metastable fee Fe has a complex antiferromagnetic state associated with a small tetragonal 
deformation. 1,4 This diversity corroborates the strong interplay between structural and magnetic degrees 
of freedom in Fe. 

The weak FM nature of a-Fe makes its magnetic interactions, and hence its characteristic properties, 
susceptible to chemical perturbations by, e.g., transition metal impurities. The effect is reflected in the 
well-known Slater-Pauling curves for the magnetic moment and Curie temperature of dilute Fe-based bi¬ 
naries. 5,6 The magneto-elastic coupling is manifested in the observed softening of the elastic parameters 
near the magnetic phase transition in bee Fe. 7 The role of magnetism in the phase stability of hydro¬ 
statically pressurized bee Fe has also been underlined. 8,9 However, explorations of the interplay among 
magnetism, lattice strain, and micro-mechanical properties are missing hitherto. 

Using computer simulations, we reveal the impact of the magneto-structural coupling on the temperature- 
dependent mechanical strength of bee Fe under large anisotropic strain. We apply uniaxial tensile load to 
a perfect bee Fe single-crystal and investigate the non-linear stress-strain response up to the point of me¬ 
chanical failure of the lattice (the maximum stress that a perfect lattice can sustain under homogeneous 
strain). The stress at the failure point under tension defines the ideal tensile strength (ITS, denoted by 
<7 m ), which is a fundamental intrinsic mechanical parameter of theoretical and practical importance. 1015 
We model the temperature effects in terms of a first-principles theoretical approach by taking into account 
contributions arising from phononic, electronic and magnetic degrees of freedom. The simulation details 
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can be found in the method section. 


Results and Discussions 

Our predicted maximum tensile strength of Fe as a function of temperature T [o m (7’)] is shown in Fig. 1 
(solid line and circles). At OK, the ITS amounts to 12.6GPa. For comparison, the measured room- 
temperature ultimate tensile strengths of bulk bcc iron ranges between 0.1 GPa and 0.3 GPa, 16 and the 
largest reported one is ~ 6 GPa for [001] oriented Fe whiskers. 17 Temperature is found to have a severe 
impact on the ITS of Fe. Namely, <7 m (T) remains nearly constant up to ~ 350K, decreases by ~ 8% 
between ~ 350 K and ~ 500 K compared to the 0 K strength and then drops by ~ 90 % between ~ 500 K 
and ~ 920 K. At temperatures above ~ 1000 K, Fe can resist a maximum tensile strength of ~ 1.0 GPa. 
For all investigated temperatures, we found that the ideal Fe lattice maintains body-centered tetragonal 
(bet, lattice parameters a and c) symmetry during the tension process and eventually fails by cleavage of 
the (001) planes. 18 

Figure 1 also shows the ITS of Fe considering merely the thermal magnetic disorder effect (dashed 
line, open diamonds). Comparing these data to the full ITS (solid line and circles), it can be inferred that 
the main trend of o m (T) is governed by the magnetic disorder term. The pronounced drop of the ITS at 
~ 500 K is related to the magnetic properties of strained Fe, more precisely to its magnetic transition tem¬ 
perature. The Curie temperature is shown as a function of the bet lattice parameters in Fig. 2. Compared 
to bcc Fe, Tq of bet Fe is strongly reduced. The main trend is that Tq decreases as the tetragonality (c/a) 
increases. A drop in Tq with increasing c/a was reported in previous fixed-volume calculations, 19 which 
turns out to be the case for the present uniaxial tensile deformations as well. 

In order to illustrate the significance of the lattice strain-induced reduction of Tq on the strength of Fe, 
we computed an auxiliary ITS [<7 1 ‘ n (7’)] assuming a constant Tq for the strained lattice. We chose without 
loss of generality the present theoretical 7 C ° of bcc Fe (1066 K) (see Table SI in the Supplementary infor¬ 
mation). The magnetic disorder effect on o^(T) is shown in Fig. 1 (dotted line, open triangles). The drop 
of the ITS o L m (T) is shifted to higher temperatures compared to <7 m (T). We conclude that the strongly 
reduced Curie temperatures of the distorted bet Fe lattices compared to bcc Fe are primarily responsible 
for the drop of the ITS at ~ 500 K. In other words, lattice deformation destabilizes the magnetic order 
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(shown schematically in insets of Fig. 2) which leads to an unexpected softening of the lattice already at 
temperatures far below . 

At temperatures T > 650 K, Fe reaches the maximum strength (fails) very close to the magnetic 
instability (T/T^ m > 0.9). The strength of these magnetically barely ordered systems approaches the 
very low strength of the PM state. The emerging question is why the ideal tensile strength of PM Fe is so 
much lower than that of FM Fe. We recall that bcc Fe is susceptible to the formation of sizable moments 
in both FM and PM states due to its peculiar electronic structure. 20,21 The nonmagnetic (NM) bcc Fe 
is thermodynamically and dynamically unstable as illustrated in Fig. 3(a), where the total energy along 
the constant-volume Bain path is shown. The pronounced peak in the NM density of states (DOS) at 
the Fermi level (£p) [Fig. 3(b)] drives the onset of FM order, that is, the majority and the minority spin 
channels populate and depopulate, respectively, causing an exchange split and moving the peak away 
from Ef. 21 When the FM equilibrium state is reached (at 0 K), £p is located near the bottom of the 
pseudo gap at approximately half band filling in the minority channel which explains the pronounced 
stability of FM bcc Fe. This is reflected by a deep minimum in the energy versus a and c map (see 
Fig. S2 in the Supplementary information). The depth of this energy minimum in connection to the ITS 
is best characterized by the total energy cost to reach the ideal strain (e m ) from the unstrained bcc phase 
(strain energy), which amounts to 0.426 mRy/atom per % strain in the FM state. 

The formation of local magnetic moments in the PM state also stabilizes the bcc structure. This is 
illustrated in Fig. 3(a), where we show the constant-volume Bain path of PM bet Fe for various fixed 
values of the local magnetic moment (/i). It is found that increasing /i gradually stabilizes the bcc 
structure and for the equilibrium value of the local magnetic moment in the PM state (2.1 /Jr), a shallow 
minimum is formed on the tetragonal energy curve. Figure 3(a) (inset) displays the total energy change of 
the bcc structure corresponding to a small volume-preserving tetragonal deformation [A£(5) — E (8) — 
£■(0)] for the same /i values as in the main figure. A E(8) is negative for NM Fe, but increasing /i turns 
A E(8) positive indicating the mechanical stabilization of PM bcc Fe upon PM moment formation. It is 
found that the mechanical stabilization due to local magnetic moment formation, embodied in the trend 
of AE(8) as a function of /i, is in fact determined by the kinetic energy (solid symbols) and thus by the 
details of the electronic structure. In the PM state, the local magnetic moment formation is due to a split- 
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band mechanism for which the average exchange field is zero . 20,2 The disorder-broadened DOSs of PM 
Fe [Fig. 3(b)] reveal that the formation of local magnetic moments effectively removes the peak at E p 
seen for NM Fe. However, this mechanism does not yield a comparably large cohesive energy increase 
as for the FM case. The strain energy cost to reach £ m starting from the shallow energy minimum equals 
0.024 mRy/atom per % strain in the PM phase, which is approximately 18 times smaller than in the FM 
case. Hence, the significantly lower ITS of PM Fe compared to that of FM Fe arises from the differences 
in the magnetism-driven stabilization mechanisms and the corresponding energy minima within the Bain 
configurational space (see Fig. S2 in the Supplementary information). 

The OK ideal strength of Fe can be sensitively altered already by dilute alloying with, e.g., V or 
Co. 2, Apart from the intrinsic chemical effect of the solute atom, there is a magnetic effect due to 
the interaction between the solute atom and the Fe host. Here we show that the strength of Fe can be 
significantly enhanced at high temperature by alloying with Co. 

We assessed the magnetic properties of the random Feo. 9 Coo.i solid solution using the same method¬ 
ology as for Fe. The alloying effect of Co on Tq is mainly connected to the strengthening of the first 
nearest neighbor exchange interactions (ii ). 5,24 We show the influence of Co on the eight first [/i( 8 )] 
and six second [/ 2 ( 6 )] nearest neighbor exchange interactions in Fig. 3(c) for the bee phase. J\ between 
both similar and dissimilar atomic species increases significantly in Feo. 9 Coo.i compared to J\ in pure 
Fe. At the same time, alloying with Co weakens Ji, but this effect is of lesser importance for Tq than the 
strengthening of J\. 

We computed the ITS of the Feo. 9 Coo.i alloy at OK 1 and in the high-temperature interval between 
500K and 1000K. As shown in Fig. 1, Feo. 9 Coo.i exhibits ~ 10% larger ITS than Fe at temperatures 
below 500 K. However, the strengthening impact of alloying on the ITS becomes more dramatic in the 
high-temperature region. The physical origin underlying this effect is the higher Curie temperature of bet 
Feo.gCoo.i compared to bet Fe, which is a consequence of both the enhanced nearest neighbor exchange 
interactions and the stronger ferromagnetism in Feo. 9 Coo.i, i.e., both the magnetic moment and exchange 
interactions are larger and more stable in the Fe-Co alloy than in pure Fe upon structural perturbation 
(tetragonalization ). 6, 23,24 The alloying effect of Co on J\ ( 8 ), .A(4), and J^{2) (the two third nearest 
neighbor interactions) is shown in Fig. 3(d) for one representative bet structure (with a — 2.732A and 
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c = 3.109 A). The addition of Co leads to a strong increase of J\ and a weakening of Ji and J 3 similar 
to the bee phase. For the chosen example, the Tcs of Fe and Feo. 9 Coo.i amount to 707 K and 1053 K, 
respectively, i.e., Co addition results in an increase of Tq which is in magnitude much more pronounced 
for the bet structure than for the bee one. 

The predicted strong magneto-mechanical softening of Fe above ~ 500 K should be observable in 
flawless systems. Indeed, for single-crystalline Fe whiskers tensioned along [001], 17,25,26 Brenner re¬ 
ported a pronounced temperature dependence of the average and maximum tensile strengths as shown 
in Fig. 1 (inset). It seems plausible to assume that the maximum strength corresponds to whiskers with 
the lowest defect density and highest surface perfection. The failure of whiskers with diameter < 6/im 
was reported to occur without appreciable plastic deformation. It is important to realize that the observed 
temperature gradient of the measured strength is much stronger than the one assuming only the thermally- 
activated nucleation of dislocations at local defects. 17 Hence, the observed softening of Fe whiskers re¬ 
quires additional intrinsic mechanisms emerging from the atomic-level interactions. The present results 
obtained for the tensile strength of ideal Fe crystal show a strong temperature gradient due to tensile-strain 
induced magnetic softening which actually nicely follows the experimental trend. The qualitative differ¬ 
ence in the theoretical and experimental temperature dependence may be ascribed to the low strain-rate 
inherent in experiment, which can allow for the activation of dislocation mechanisms. 27,28 

Conclusions 

We have discovered a strong magneto-mechanical softening of iron single-crystals upon tensile loading. 
We showed that the strength along the [001] direction persists up to ~ 500 K, however, diminishes most 
strongly in the temperature interval ~ 500 - 900 K due to the loss of the net magnetization upon uniaxial 
strain. The strength in the paramagnetic phase is more than 10 times lower than that in the ferromagnetic 
phase. We found that Fe fails by cleavage at all investigated temperatures. 

We have demonstrated that the intrinsic strength of Fe at high temperatures is significantly enhanced 
by alloying with Co. This finding opens a way to carefully scrutinize the proposed connection 29,30 be¬ 
tween the tensile strength of a defect-free ideal crystal and the measured tensile strength of single-crystal 
whiskers or nanoscale systems (e.g., nanopillars). Measuring the tensile strength of Fe-Co whiskers at 
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elevated temperature could verify the here predicted differences in the high-temperature intrinsic mechan¬ 
ical properties of Fe and Fe-Co alloy. 

Extending the above mechanism to other systems, we expect that the large magnetic effect on the 
strength exists for all dilute Fe alloys and that the magneto-mechanical softening temperature can be 
sensitively tuned by proper selection of the alloying elements. Solutes that strengthen (weaken) the 
ferromagnetic order in the tetragonally distorted Fe increase (decrease) the intrinsic strength of defect- 
free Fe crystals at temperatures above (below) ~ 600 K. 

Method 

The adopted first-principles method is based on density-functional theory as implemented in the exact 
muffin-tin orbitals (EMTO) method 31-33 with exchange-correlation parameterized by Perdew, Burke, 
and Ernzerhof (PBE). 34,35 EMTO features the coherent-potential approximation (CPA), which allows 
to describe the disordered PM state by the disordered local moment (DLM) approach. 31 

The ITS (<7 m ) with corresponding strain (e m ) is the first maximum of the stress-strain curve, o(e) = 
VJFj^ ( e ) an d C(e) are the relaxed volume and the calculated free energy at strain e], upon uniaxial 
loading. The tensile stress was determined by incrementally straining the crystal along the [001] direction 
(the weakest direction for bee crystals) and taking the derivative of the free energy [G(e)] with respect 
to e. The two lattice vectors perpendicular to the [001] direction were relaxed at each value of the 
strain allowing for a possible symmetry lowering deformation with respect to the initial body-centered 
tetragonal symmetry. 3 - 1 

The temperature induced contribution from electronic excitations to the ITS was considered by smear¬ 
ing the density of states with the Fermi-Dirac distribution. 31 To measure the effect of explicit lattice 
vibrations, we computed the vibrational free energy as a function of strain in the vicinity of the stress 
maximum within the Debye model 7 employing an effective Debye temperature which was determined 
from bulk parameters. 38 

In order to take into account the effect of thermal expansion, the ground-state bee lattice parameter 
at temperature T was obtained by rescaling the theoretical equilibrium lattice parameter with the experi¬ 
mentally determined thermal expansion. 39,40 The expanded volume was stabilized by a hydrostatic stress 
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Pt derived from the partial derivative of the bee total energy with respect to the volume taken at the ex¬ 
panded lattice parameter corresponding to T. pr was accounted for in the relaxation during straining the 
lattice. To this end, the relaxed lattice maintained an isotropic normal stress pr in the (001) plane for 
each value of the strain. 

We modeled the effect of thermal magnetic disorder on the total energy by means of the partially 
disordered local moment (PDLM) approximation. 41,42 Within PDLM, the magnetic state of Fe and 
Feo. 9 Coo.i alloy are described as a binary Fc{_ r Fei and a quaternary (Fe/, 9 Coq | )i_ x (FeQ 9 Co ( j { ) x alloy 
with concentration x varying from 0 to 0.5 and anti-parallel spin orientation of the two alloy components. 
Case x — 0 corresponds to the completely ordered FM state with magnetization m = 1. As x is gradually 
increased to 0.5, the magnetically random PM state lacking magnetic long and short range order (m = 0) 
is obtained (DLM state). The PDLM approach describes the energetics underlying the loss of magnetic 
order. Connection to the temperature is provided through the magnetization curve [m(z), z = T/Tq 
being the reduced temperature] which maps the computed m to T. It is important to realize that for Fe 
and Feo. 9 Coo.i the thermal spin dynamics embodied in the shape of m(z) and in the Tq value change 
under the presently applied uniaxial loading. The computational details for Tq and m(z) are given in the 
Supplementary information. 
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Temperature (K) 

Figure 1 . The ITS of bcc Fe in tension along the [001] direction as a function of temperature (<7 m ('/’), 
solid line and circles); the ITS taking into account only magnetic disorder (dashed line, open diamonds); 
the ITS taking into account magnetic disorder and neglecting the change of Tq with structural 
deformation (o^, dotted line, open triangles). The Fco.gCoo, i alloy (stars) possesses a slightly larger ITS 
at 0 K, but compared to pure Fe the ITS drastically increases at high temperatures. All data are 
normalized to the ITS of Fe at 0 K (12.6 GPa). (Inset) Experimentally determined tensile strength of 
[001] oriented Fe whiskers as a function of temperature from Ref. 1 The two data sets give the average 
tensile strength and the maximum tensile strength for whiskers possessing nearly equal diameter 
(5.1-5.4/im). Lines guide the eye. 
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Figure 2. Contour plot of Tq of bet Fe as a function of the lattice parameters a and c. The region where 
7c is shown is confined by the bcc ground state (black asterisk) and the failure points (e m ) at different T 
(stars, numbers denoting T). The Curie temperature corresponding to each failure point (7 )?'") can be 
read from the legend. The ground state bcc structure possesses the highest calculated Tq in the region 
shown. The dashed line represents the hyperbola of constant volume equal to the bcc equilibrium 
volume. The insets sketch local spins for the unstrained parent lattice and for the strained lattice at high 
temperature, illustrating the magnetic disorder increase with tensile strain. 
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Figure 3. (a) Total energy of NM bet Fe and PM bet Fe for various values of the local magnetic 
moment. The energies are plotted with respect to the energy of the PM bcc state with equilibrium local 
magnetic moment (2.1 /1 b), and the volume is fixed to that of PM bcc Fe. The inset displays the total 
energy change (open symbols) and the kinetic energy (solid symbols) change corresponding to a small 
volume-preserving tetragonal shear (arbitrary units), (b) Total DOS of PM bcc Fe for different local 
magnetic moments, (c) and (d) The exchange interactions J\ and Ji of bcc Fe and Feo. 9 Coo.i and J\-J 3 
of one representative bet structure {a = 2.732 A, c — 3.109 A). 
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